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Abstract 

The auxiliary linear problems are presented for all discretization 
levels of the Hirota-Ohta system. The structure of these linear prob- 
lems coincides essentially with the structure of Nonlinear Schrodinger 
hierarchy. The squared eigenfunction constraints are found which re- 
late Hirota-Ohta and Kulish-Sklyanin vectorial NLS hierarchies. 
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1 Introduction 

,__! The Hirota-Ohta system [13] was introduced by means of Pfamanization 

applied to the Kadomtsev-Petviashvili equation, that is as a result of a cer- 
tain procedure which allows to replace multi-soliton solutions represented 
by determinants with solutions represented by Pfaffians. Later on, this pro- 
cedure was applied to many other equations, see e.g. [10]. The differential- 

t — difference member of the Hirota-Ohta hierarchy known as the Pfaff lattice 

[1, 2, 3, 16] was introduced within the theory of random matrix models. 
The whole hierarchy can be derived also within the general approach based 
on Clifford algebra representations and the boson-fermion correspondence 

•»h [14, 15], for example, this technique was used for the derivation of Backlund 

rN transformation [201. 

In this Letter we suggest an alternative and more simple derivation of 
the Hirota-Ohta hierarchy by use of Zakharov-Shabat method [26] based on 
the study of auxiliary linear problems. To the best of our knowledge, such 
linear problems for the continuous part of the hierarchy appeared so far only 
in paper [17]. The completely discrete Hirota-Ohta system was found for 
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the first time, apparently, in paper [12]. However, it is not clear from this 
work, whether this discrete system belongs to the Hirota-Ohta hierarchy or 
just serves as its discrete analog. In Section 3, we answer in the affirmative 
to this question and provide the consistent linear problems for all discrete 
levels of the Hirota-Ohta system, that is for the systems with one, two and 
three discrete variables. 

The linear problems under consideration are rather deep reductions. For 
example, in the completely discrete case, we consider a 2-component 4-point 
scheme with three field variables only, while the generic scheme of such type 
corresponding to the discrete Darboux-Zakharov-Manakov system [22, 7] 
contains eight fields (two 2x2 matrices). The following empirical observation 
allows to recognize these reductions: 

the structure of linear problems for 3-dimensional Hirota-Ohta 
hierarchy coincides with the structure of 2-dimensional hierarchy 
of Nonlinear Schrodinger equation (NLS). 

The difference is that the nonlinear terms are replaced with the linear ones 
with variable coefficients, the Hirota-Ohta hierarchy arising from the com- 
patibility conditions for these coefficients while the NLS flows commute iden- 
tically. 

In Section 2 we recall some standard equations from the NLS hierar- 
chy: the third order higher symmetry, Backlund transformations and the 
nonlinear superposition principle. These equations are used as a hint while 
choosing the form of the linear problems for the Hirota-Ohta hierarchy in 
Section 3. A more precise conformity with the theory of 2-dimensional sys- 
tems is pointed out in Section 4 devoted to the vectorial generalization of 
NLS by Kulish and Skkyanin [19]. This correspondence was found in [5] by 
one of the authors unaware of Hirota-Ohta results at that moment. 

2 NLS hierarchy 

The following notations are used below. The field variables ip, <fi depend on 
the infinite set of independent continuous variables t\ = x, t-z = y, £3 = t, . . . 

and integer- valued variables m, . . . , rii, . . . , rij, . . . , i%k, Subscripts x, y, t 

denote the partial derivatives, and subscripts l,i,j,k denote the shifts on 
the discrete variables: ipi = ip(. . . ,n{ + 1, . . . ), minus sign denoting the 
backward shift: ip-i = ip(. . . , rii — 1, . . . ). 

The Nonlinear Schrodinger equation is of the form 

i> y = tpxx + 2V> 2 </>, -4> y = (p xx + 24> 2 ip. (1) 

Its simplest higher symmetry is of third order: 

l\>t = '^xxx + ( o4>^x- l (f>t = <l>xxx + 6i><i>(f>x- (2) 



The discrete variable n\ will play the distinguished role in what follows. 
The corresponding shift defines Backlund-Schlesinger transformation as the 
explicit mapping 

•>h = Ax - 4> 2 J4> + i> 2 fa <h = i/^ (3) 

which acts on the solutions of systems (1), (2). The iterations of this map- 
ping are governed, under the change ip = e q , <f> = e -9-1 , by the Toda lattice 

qxx = ^- q -e q - q - 1 - (4) 

All other discrete variables correspond to Backlund-Darboux transforma- 
tions; x-parts of these transformations are of the form 

t/j x = <;/>; + a (< V + ^ 2 fa, -<j> i , x = (/) + a ii) <t>i + (l> 2 ip, (5) 

where or- 1 ' is an arbitrary parameter, associated with i-th discrete coordinate 
direction and depending on ni only (that is, cr- = «w at j ^ i). Second 
equation (5) defines fa as a solution of Riccati equation, then fy is explicitly 
found from the first equation. A generic Backlund transformation for the 
NLS equation is decomposed as a sequence of elementary transformations 
of the form (3), (5) and their inverses. 

The completely discrete part of the NLS hierarchy is described by 5-point 
equations of discrete Toda type 

e ii,-i-i _ e 9-q-i,i + e ii-i _ e 9-q-i + a W _ a W = o (6) 

and 2-component quad-equations 

^j = ^+ K — 7— — , fa = fa+ y —. tt^- 1 - 7) 

1 - Ipfaj 1 - ^0jj 

Equations (6) define the permutability property of the transformations (3) 
and (5). This means that equation (6) is consistent with the differentiation 

q x = e *-9 + e ?-9-i,« + a W = e ft,-i-9 + e ?-«-i + a^, 

that is the derivative with respect to x of the left-hand side of (6) vanishes 
in virtue of (6) itself. Moreover, the variables q(n\,rii) at fixed n, satisfy 
the Toda lattice (4), and variables ip = e q ( ni,ni >, 4> = e~ q ( ni ~ 1,ni > at fixed n\ 
satisfy the lattice (5). 

Analogously, equations (7) define the permutability property of two 
transformations of the form (5). Let \& = (ip, fa), then the following property 
holds: if \P and \Pj are related by the transformation (5) with parameter ay\ 
and ^i and ^ij are related by the transformation with parameter a^' then 
the variable \Pj defined accordingly to (7) is as well related to ^ and Vl/y by 





Figure 1: RjkRikRij — RijRikRjk 

the transformations of the same type, but with the interchanged parameters 
cr J ) and aW. 

The mappings Rij : (fy, \Pj, ^ij) *- > ^j satisfy an important 3D-consistency 
property expressed by the Yang-Baxter identity (see fig. 1) 



^jk^ik^ij — ttijrCik-tijk. 



(8) 



An essential distinction with the scalar quad-equations [6] is that equations 
(7) cannot be solved with respect to variables ^f or ^fij. Hence the formu- 
lation of 3D-consistency is possible here only under a choice of the initial 
data along some sequence of the vertices of the cube, for example *&, ^i, 

^12, #123- 



3 Hirota-Ohta hierarchy 

We keep the previous notations for the derivatives and shifts with respect 
to the independent variables. As for the dependent variables, now ip and 
4> play the role of wave-functions, while the field variables will be denoted 
u, v, w, u)W , w^\ q (coefficients of linear problems) and /, g, h (r-functions). 
Superscripts in w^ l \ w^' indicate that these variables are associated with 
the corresponding direction or plane in the lattice (but, in contrast to pa- 
rameters a^' from the previous section, these depend on the whole set of 
independent variables). 

We accept the following linearization of the NLS equation (1) as a start- 
ing point: 



Ipy = Ipxx +Wlp + 2u(j), -(j)y = 4> xx + 2vip + W(f). 

Now, let us search for a commuting flow of the form 

ipt = ipxxx + aipx + bip + ccf), 4>t = 4>xxx + Mx + Bip + C(j) 



(9) 



with undetermined coefficients (that is, the cubic terms in (2) are replaced 
with the linear combination of the factors). An easy computation specifies 
the coefficients: 

3 3 

fpt = ^xxx + ^wtp x + -(w x + q)ip + 3u x cf), 

3 3 

4>t = 4>xxx + -jW&x + j{w x - q)4> + 3v x ip (10) 

and yields the compatibility conditions of equations (9) and (10) in the form 
of the following system: 

— 2ut = u xxx — 3u xy + 3wu x — 3qu, 
-2v t = v xxx + 3v xy + 3wv x + 3qv, 
4w t = w xxx -24(uv) x + 6ww x + 3q y , q x = w y . (11) 

This is nothing but the Hirota-Ohta system [13]. It is clear that its higher 
symmetries can be derived quite analogously by use of the linear problems 
of higher order with respect to d x ; however, we will skip this. 
The change of variables 

n=-, v = -, w = 2(logh) xx , q = 2(logh) xy (12) 

brings the system (11) to the bilinear form 

(2A - 3D x D y + D'Df -h = 0, 
(2D t + 3D x D y + Dl)g -h = 0, 
{AD x D t - 3D 2 y - Di)h ■ h + 24/ 5 = (13) 

where D z a- b = a z b — ab z denotes the Hirota operator. The variables f,g,h 
are especially convenient for rewriting of the discrete dynamics. 

It is natural to search for an analog of the transformation (3) in the form 

V>i = aijj xx + bij) x + ap + d(f>, 0i = pip. 

The answer is not unique, because of the gauge transformations 

tl> = an(y)i>, (j)=——<i), u = — jU, v = —-v, w = w 

which leave invariant the form of basic linear problem (9). It is possible 
to fix some coefficients by use of these transformations and to choose the 
mapping in the following form: 

i>i = i>xx -ipx + ( w + -^- y -)ip + u(j), <f>i = V- (14) 

u V 2u / u 

Its action on the coefficients of the linear problems is described as follows. 



Statement 1. System (11) admits the explicit auto-transformation 

2 / \ ^ / 2\ / 2 

Ul = U V + -(U X W X - UWy) -\ [UU XX - U x ) + —{UUyy ~ U y 

o \ 

ZUU XX y -\- ZU x ll X y -\- uu xxxx zu x u xxx -\- U xx j , 
vi = l/u, wi=w + 2(logu) xx , qi=q + 2(\ogu) xy . (15) 

The substitutions (12) reduce the last three equations in (15) to h\ = f, 
gi = h, that is the iterations of this mapping generate the sequence 

... / = /i(ni + l), h = h(n 1 ), g = h{m - 1) ... 

Moreover, the first equation (15) is equivalent to the relation 

A \ Bi B 
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.h\h' x h\ h 2 ' 

where 

A = (3D x D y - D x )hi -h, B = {3D 2 y + D x )h ■ h - 24hih-i. 

This relation is not bilinear, however it turns into identity in virtue of the 
system (13) which takes the form (cf [1, 2, 3]) 

(2A - 3D x D y + D*)hi -h = 0, 

{AD x D t - 2,D 2 y - D x )h ■ h + 24hih-i = 0. (16) 

Indeed, in virtue of these equations 

A B h t 

MT 2ri - 2r ' ¥ = 8r - r = J- 

Darboux-Backlund transformation is derived from the linear problem 

ipx = ipi + W^1p + U(l>i, -<$>i,x = 4> + Vilj) + W^(j)i, (17) 

and an analog of the nonlinear superposition principle (7) is described by 
the linear problem 

i/fj = tpi + w {ij) {^ + uuffij), 4>j = &- w^Xviji}) + <pij), i < j. (18) 

The compatibility conditions for these equations define all discrete levels of 
the Hirota-Ohta system. We write them in the bilinear form at once, since 
the equations for the coefficients are rather bulky. To do this we use the 
following relations, in addition to the substitutions (12): 

W » = ^-%, w^ = ^. (19) 

hi h hihj 



Statement 2. The compatibility conditions for equations (9) and (17) yield 
the system 

hgi, y - h y gi + hg iyXX - 2h x g i<x + h xx gi - 2hig = 0, 

hh i:V - hyhi - hh i)XX + 2h x h itX - h xx hi + 2/gj = 0; (20) 

the compatibility conditions of (18) and two equations of the form (17), 
corresponding to i and j, bring to the system 

JilijjX Jx'l'ij — Jiilj Jji^ii 
^gij,x ft-xgij — <H9j FljQii 

hih j)X -h i>x hj = fgij - hhij, i < j; (21) 

the compatibility conditions for three linear problems of the form (18), cor- 
responding to i,j and k, yield the system 

J hijk Ji<ljk t Jjh"ik Jkfcij = U, 

hgtjk — higjk + hjgik — h k gij = 0, 

fgijk - Khjk + hjh ik - h k h i:j = 0, i < j < k. (22) 

Equations (22) are nothing but the discrete version of the Hirota-Ohta 
system found in [12]. Let us explain its derivation in more details. Formally 
one may consider the linear problem of form (18) as a very special case 
of 4-point scheme for 2-component function \P = (-*/>, </>). As in the case of 
the NLS equation, <j) and ipy are not involved in the equations, and this 
dictates the following choice of the initial data needed for the computation 
of compatibility conditions: \1/, \I/j, *$>ij, *&ij k . For the sake of definiteness we 
will assume that i < j < k. Then the compatibility condition is expressed 
by Yang-Baxter equation (8) where Rij is the mapping (\£, ^j, ^>ij) h- >■ 'Un- 
defined by equations (18). The computations use the following equations: 

ipj = ^i + w^ (V> + u<t>ij), ifjjk = ipik + Wfc (V'fe + u k 4>ijk), 



R 



V ' rh, = rh, - 7/)(«)Ci),,7/) 4- A.A A J... „.,(*i) 



W 



{Vijlp + 4>ij ) , (f> jk = (f) ik - w ( v> (v ijk 1p k + (j> ijk ), 



i>jk = ipij + w? (i>j + Ujfajk), ifj k = ipi + w^ (ip + ucj) ik ), 

{tk \vijkipj + <f>ijk), 4>k = $i- w^iviki; + 4> ik ) 



rLi k . i, 

4>jk = 4>ij - w) ' [Vijk'tPj t- 9ij 

. Ipk = Ipj + W^ k ~> (ip + U(f>j k ), Ak = i>ij + Wl jk > {ipi + Ui<j) ijk ), 

1 '" v ' (j) k = <f>j - w^ (v jk iP + 4> jk ), (j) ik = fcj - W W ( Vijk ^i + <t> ijh ). 

Left hand side of Yang-Baxter equation (8) corresponds to the consecutive 
computations according to the left column downwards, while the right hand 
side corresponds to the right column computed upwards. The resulting non- 
linear relations for the coefficients contain, in particular, the conservation 
laws which imply the parametrization of w^> given in (19). The rest rela- 
tions then turn out to be equivalent to system (22). 



4 Kulish-Sklyanin vector hierarchy 

A wide and well-known class of reductions from 3-dimensional equations 
to vectorial 2- dimensional ones consists of so-called squared eigenfunction 
constraints. For instance, the Manakov system [21, 9] and its third order 
symmetry 



ipy = ^xx + 2(ip, 4>)i>, -<j) y = (f) xx + 2(ijj, 
l/>t = Ipxxx + 3(V>, 4>)^x + 3{lpx, (f>)4>, <t>t = 4>xxx + 3(ip, <P)4>x + 3(V>, (f>x)4> 
define such a reduction for the Kadomtsev-Petviashvili equation 
Aut = u xxx - 6uu x + 3q y , q x = u y 

with respect to the quantities u = —2{tp, tfi), q = 2{ip, 4> x ) — 2{ifi x , 4>) [18]. 

Now we demonstrate that the Hirota-Ohta hierarchy (for all discrete 
levels) arises analogously from another vectorial version of the NLS equation 
introduced by Kulish and Sklyanin [19]: 

Ipy = Ipxx + 4(V>, (j))lp ~ 2(-0, 1p)4>, 

-0„ = ^ + 4(^,^-2(0,0)^. (23) 

Let us write the corresponding generalizations for the basic equations in 
Section 2. The third order symmetry reads 

4>t = Ipxxx + 6(V>, (P)lpx + 6{lpx, 4>}lp - 6 (V ; , 1px}(f>, 

<h = <f>xxx + 6{ip, 4>)4>x + 6(i/), 4>x)4> - 6(0, (j) x )ip- (24) 

The analog of Backlund-Schlesinger transformation (3) is of the form [25] 

1* = A* - * { 4r£** + { 4r^r^ + 2(^, M - (A ^ 

V, (25) 



the analog of Backlund-Darboux transformation (5) reads 

4>x = i/>i + t* (i ty + 2(V>, <k)ij> - {il>, ip)<f>i, 
-<Pi,x = <P + a^<j)i + 2(ij>, <j>i)<t>i - (fa, <tn)ij) (26) 

and, finally, the analog of nonlinear superposition principle (7) is [4] 

, = , _ («»-q^)(^--(^-,^W ( 7) 

The comparison with the linear problems from the previous section proves 
immediately that the Kulish-Sklyanin hierarchy is their self-consistent re- 
duction. More precisely, the following statement holds true. 



Statement 3. Equations (23)-(27) are consistent and, in virtue of these 
equations, the quantities 






,w) _ ° w - ggj 



1 - 2(ip, fcj) + (V>, ip){(f)ij,<j)ij) 
satisfy the equations of the Hirota-Ohta hierarchy. 

5 Concluding remarks 

The main result of this Letter is an uniform derivation of continuous and 
discrete equations of the Hirota-Ohta hierarchy from the compatibility con- 
ditions of auxiliary linear problems. The structure of these problems pat- 
terns after the structure of 2-dimensional NLS hierarchy. We hope that this 
observation may be useful for other 3-dimensional equations as well. 

It is worth noticing also the relation of above discrete linear problems 
with self-adjoint 5- and 7-point schemes, on square and triangular lattices 
respectively. The linear problems of this type attracted much attention 
recently. In particular, it is known that they appear from the scalar 4- 
point schemes on square and rhombic lattices as a result of restriction on 
even/odd sublattices [8]. In our case, such linear problems appear from 2- 
component linear problems presented in Section 3 as a result of elimination 
of one component. For instance, the elimination of <fi from (14) and (17) 
yields the following 5-point equation (compare with discrete Toda lattice 
(6)): 

-fa + — i/>-i - -ih,-t - — V>-m + - W°_i + ™ {i) - — ) V = o. 

U U-i u M— i,i u\ ' u J 

Analogously, eliminating <p from equations (18) results in the 7-point linear 
problem in the plane (ij). We mention that some examples of continuous 
Toda-like dynamics associated with 5- and 7-point linear problems were 
studied in the papers [24, 11, 23]. It remains unclear, whether these examples 
coincide with the flow (17) of the Hirota-Ohta hierarchy; more likely, they 
may interpreted as its negative flows. 
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